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SUMMARY 


Pratt & Whitney Aircraft Group, Government Products Division (P&WA/GPD) con- 
ducted this program under National Aeronautics and Space Administration (NASA) Contract 
NAS;i 21607 to investigate techniques for improving turbofan engine digital model 
high frequency response A previous contract. NASH-20292, optimized the dynamic configura- 
tion ot a TFBO engine model using strictly one-dimensional modeling techniques with limited 
success in matching the response of the engine data. This program developed a fan model 
using i wo dimensional modeling technique* to add dynamic radial coupling between the core 
stream and the bypass stream of the fan. When incorporated into a complete TF30 engine 
simulation, the fan model greatly improved compression system frequency response to planar 
inlet pressure disturbances up to 100'Hz. The improved simulation also matched engine 
stability limits at l.THz, whereas the one-dimensional fan mi>del required twice the inlet 
pressure amplitude to stall the simulation. With verification of the two-dimensional fan model, 
this program formulated a high frequency FlOOtH) engine simulation using row-by-row com- 
pression system characteristics. In addition to the FI (HUH) remote splitter fan, the program 
modified the model fan characteristics to simulate a proximate splitter version of the FlOO(H) 
engine. The response of the FI (Hi models was compared to high frequency (50-60 Hz) 
augmeiitor rumble data with good comparison of phase, but differences in gain. The source of 
tlie FIDO model inaccuracy was traced to the fan duct model. 


INTRODUCTION 


(Jas turbine engine system design stresses the importance of developing the capability to 
model high frequency phenomena and, thus, properly simulate the effects of augmentor and 
inlet disturbance* on compression system response and stability. For example, augmentor 
transient disturbances, which have been tracked forward through the fan duct and around the 
inlet splitter, produce destabilizing effects on the compression system. Design studies indicate 
that transient operation of fan duct burners, variable geometry, and air invertor valves also 
result in duct flow disturbances that impact compression system stability. Evaluation of 
dynamic engine stability early in the design process requires a good dynamic engine model to 
ensure propulsion system operational suitability. 

NASA C'ontract NASH- 20292 developed a high-frequency TF30 turbofan engine model 
with the objective of achieving reponse characteristics up to 100Hz. Modifications to the 
compression system improved the dynamic response characteristics of the engine simulation 
based on an investigation of the influence of size and number c f finite elements using strictly 
one dimensional modeling techniques. Based on the results of this program, high-frequency 
engine simulations using finite elements can be assembled to optimize the engine dynamic 
configuration with respect to dynamic characteristics and computer execution time. Resizing 
the compression system finite elements improved the dynamic characteristics of the engine 
simulation, hut showed that additional refinements were required to obtain close agreement 
between simulation and actual engine dynamic characteristics. 

Modeling of the fan bypass/core stream dynamic coupling and the transient splitter 
aerodynamics proved to he the main areas of needed improvement. Previous TF.'M) model 
results indicated a pressure disturbance must be twice the amplitude observed in tests to 
generate model instability. The simulation showed the fan to be the stability-limiting compo- 
nent where test data shows rotating stall to he initiated in the low-pressure compressor. The 
frequency response at the fan discharge of the model did not match recorded engine data 
above approximately 20*Hz. The gain of fan discharge pressures in response to spatially 
uniform (planar) periodic pressure disturbances introduced at the engine inlet did not compare 
as well as desired with engine data. 

The technical approach of the turbofan model frequency response improvement program 
conducted under contract NAS3 21(507 sought to provide dynamic radial coupling in the fan by 
describing the gas path flowfield of the fan in a two-dimensional format of radial as well as 
axial directions. Three describing equations of continuity, energy, and momentum across a 
finite element were provided in the radial direction between the finite elements representing 
the fan core stream and the fan bypass stream. 

This two-dimensional fan concept was constructed into a computer model and integrated 
into a complete turbofan engine model. Using TF30 engine component characteristics, the 
frequency response of the turbofan model was compared to measured engine data. When 
significant frequency response improvements resulted in the new TF30 model compared to the 
one-dimensional model, the FIOOCO engine was then modeled. In addition to the baseline 
remote splitter FIDO model, a proximate splitter Fl(M) model was also constructed. The 
response of both F100 models was compared to measured engine test data. 


TWO-DIMENSIONAL PAN MODEL 


UajIaI 

MODVI vonofpi 

A majority of turbofan operation occur* with a radial warpage in discharge pressure 
arrow the span of the fan. The sensitivity of the turbofan cycle balance requires a fan model to 
account for this pressure warpage during all operating conditions. An average bypass fan 
section and an average core fan section, as shown in figure 1, accomplishes this in the 
cost effective, one dimensional gas path modeling technique. A great deal of accuracy has been 
obtained in steady-state simulations by proper correlation of this warpage with some parame- 
ter. such as bypass ratio. 
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However, steady state correlations of warpage. which couple the bypass fan representa- 
tion with the core fan representation, are not necessarily sufficient for a dynamic represents 
tion Determination of transient response for stability and suitability investigations requires 
definition of accurate transient flow paths for components in the model. Proper definition of 
these paths, symbolized bv the shift in bypass and core dividing streamlines as a function of 
time, further requires the incorporation -if radial dynamic coupling in the fan component. 

This refined modeling technique allows radial communication between control volumes 
that describe the inside diameter (ID) and outside diameter (OI)) characteristics of the fan. In 
a previous NASA contract (NASH 2025)2). the dividing streamline, based on the steadvstate 
operating bypass ratio of the fan. set the characteristics of the control volumes for the bypass 
and core streams. However, no direct communication path existed between the control volumes 
in the radial direction during a transient. The new model accomplishes this communication by 
solving the basic unsteady flow equations for the radial direction in additional to axial 
direction calculations, as shown in figure 2. 
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Row ChorodoritM©* 


The one dimensional fan model utilized th* row-by-row charactariatics of NAS3- 1863ft aa 
a hflKe for the TF30 engine and NAS't 310 aa a baa* for th* F100 angina. Tha** empirical 
characteri»»'t> describe the fan performance in th* core (ID) atraam and in the bypaaa (OD) 
stream. A ..viding" streamline, aa shown in figure 3, aaparataa tha core stream and bypaaa 
stream Since the flow par unit area remain* approximately constant separate OD and II) 
pressure rise characteristics were described for each rotor and stator, independent of the 
streamline position These characteristics, shown in Appendix B for the F100 fan. utilized the 
following definitions: 
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The total and static pressures required for the calculations of 4> and 'F are related through 
Mach Number 
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where: 

A flow area 

M Mach number 

N rotor speed 

1* pressure 

T temperature 

A pressure, temperature differential 

h pressure standard pressure 

ii temperature/standard temperature 


Subscripts: 


'I' 


static 

total 





Two-Dimtnsional Equation* 


The mathematical abroach to the two-dimensional flowfield follows as a natural 
extension of the basic one-dimensional flow through a finite element, as described in Appendix 
A. Conceptually, the elements parallel each other with a primary flowpath in the axial 
direction and a secondary flowpath in the radial direction. This generates the grid that makes 
up the boundaries of the finite elements, as shown in figure 4. The velocity at the entrance and 
exit of the element contains both axial and radial components. A four-point averaging of two 
adjacent elements enables the computation of velocity at the radial boundary between the 
elements. The basic flow equations required to define the unsteady flow in an element are 
continuity, momentum, ana energy. Each equation describes a basic property of the fluid flow: 
namely mass, momentum, and energy. It should be noted that the properties of mass and 
energy within the finite element are scalar quantities, while momentum is a vector quantity. 

The following equations for time derivatives of the properties can be written for Element 
A in figure 4. 

Continuity Equation 
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Where: 

V axial velocity component 

V radial velocity component 


Momentum Equation 
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Knotty Equation 
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Subscripts 

A 

r 

Tui 
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i 1 % i-.s >n t C »liul vj im 

specific In* it at constant pressure (function of T) 

specific heat at constant volume 

energy in the element 

tone »»n t he thud due to turbomachinery 

nuns m lb; element 

moment inn in tin element 

statu pressure 

heat addition across the element boundaries 

term to an-oimt tor turhomaelnnery radial wnrpage 

radius 

statu* temperature 
voUi.it \ 

mechanical \v.»rk 
density 


element A of figure 1 
radial direction 

turbomachinery tin its inside element 
axial direction 

nodal coordinates ot figure 1 
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The continuity equation define* the time rate of change of maaa in the element in terms 
of the flow m and out of the element. The momentum equation expresses the time rate of 
chuuge of momentum in the element due to mass transfer and forces on the fluid in the axial 
direction, us well as the radial direction of flow. The energy equation represents the time rate 
«l change of energy as the sum of the kinetic and potential energies of the fluid, the rate of 
mechanical work done on the fluid, and the rate of heat addition to the fluid. Expressed in an 
integral form, the equations define the quantities of mass, momentum, and energy within the 
finite element in the following manner: 

I'nntmuitx Equation 

M { /.dt Volume) ) />(x)A(: • • (11) 

The equation integrates by assuming a linear distribution of properties across the finite 
element as follows; 

iiiii) • 1/1 * dx (12) 
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Thus, the continuity equation becomes 
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which sutifios the continuity principle; 
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Similar expressions have been derived for (he momentum and energy equations. Equa- 
tions 1 1 through 21 define the fluid conditions in Element A of figure 4. 
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Radial Term 

Because the one dimensional core stream and bypass stream eharaet eristics represent 
average pressures in their respective streams, as shown in figure 1, a radial profile in pressure 
may exist at a steady state operating point. An increment in radial momentum flux (K* in the 
above equations) balances the effect of this pressure profile. Kffort was extended without 
success during the contract period to analytically predict this radial term as a function of 
blading design and operating conditions. Thus, the model results reported herein consist of 
calculating the required steady state radial momentum term and retaining the vector angle 
between the radial and axial momentum fluxes as a constant during transients. For small 
transients away from a steady state point, this method is probably adequate. However, a new 
technique should he developed for large excursions. 
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TF30 ENGINE MODEL 


The TF30 is a twin -spool augmented turbofan engine that mixew bypass and core gas 
flows in the afterburner duct. The low-pressure spool is composed of a three-stage fan and a 
six-stage compressor driven by a three-stage turbine. The high-pressure spool consists of a 
seven-stage compressor driven by a one-stage, air-cooled turbine. 1 * combustor is a can type 
surrounded by a secondary chamber that uses high-pressure compressor discharge air to cool 
the burner cans and dilute the fuel-air mixture. Surge bleeds are incorporated at the seventh 
and twelfth stage of the compression system 

Engles Teat Data 

NASA Lewis Research Center (LeRC) tested a TF30-P-3 engine to investigate engine 
dynamics. These high frequency response results enabled an evaluation of the dynamic 
characteristics of the engine simulation in the contract. Testing occurred in an altitude 
chamber using the LeRC air jet system described in Reference 2, which induced dynamic 
pressure variations at the front end of the engine. The average inlet conditions for the test 
produced a Reynolds number index of 0.5 at an inlet pressure of 50.7 to 54.5 kPa abs (7.35 to 
7.90 psia). Kngine exhaust pressure was set at 16.5 kPa abs (2.4 psia) for all tests. As 
illustrated in figure 5, high response total presson* instrumentation locations included the 
engine air inlet tSta. 2i0), the fan first stator passage (Sta. 2.1), the fan core stream discharge 
(Sta. 2.3). fan bypass stream discharge (Sta. 2.3F), low-pressure compressor discharge (Sta. 
3.0), twelfth stage stator passage (Sta. 3.12), and high pressure compressor discharge (Sta. 4.0). 

NASA devoted a portion of the test program to an investigation of the engine dynamic 
response to spatially uniform (planar) periodic pressure disturbances induced at the engine 
inlet. By using the frequency sweep test technique described in Reference 2, gain and phase 
characteristics of response to the engine inlet pressure were established for pressures at the 
internal engine stations of 2.3, 2.3F. 3.0, and 4.0. Figures 6 through 9 akow the Sta. 2.3 and Sta. 
2.3F gain and phase characteristics of the engine operating at 8600 Nl/\fl used for com- 
parisons in this contract. Although the air jet stream was mechanically capable of frequencies 
up to 200 Hz, Reference 2 cautions against using the data above 80 Hz because of the 
measurable amplitudes in the engine inlet and fan-compressor system pressures. Investigations 
also included engine stability limits during the NASA test program. During these tests, 
discrete frequency inlet pressures were induced and amplitude increased until either 
stall/surge was encountered or the amplitude limits of the air jet system were reached. Results 
of these tests in terms of P amplitude at stall for 15 Hz frequency, as shown in figure 10, 
provided data for comparison to the stability limits of the simulation. 

Batalin* Element Configuration 

This modeling technology improvement program investigated methods to improve tur- 
bofan engine system model response in the frequency range up to 1(H) Hz. A previous NASA 
contract (NAS3-20292) utilized the TF30 mixed flow turbofan engine in strictly 
one-dimensional simulation technology studies. Through reduced-order analysis of the com- 
pression system simulation element breakup, NAS3-20292 showed that the dynamic character- 
istics of a stage-by-stage simulation at ltK) Hz could be obtained with a more simplified 
dynamic configuration, while realizing substantial savings in computer execution time. The 
analysis showed that the dynamic elements should be uniformly sized to achieve the best 
frequency response. The optimum dynamic element configuration for the TF30 engine under 
this program was 27 elements. This became the baseline for the present contract to investigate 
simulation improvement by utilizing two-dimensional effects in the fan. 
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The hemline 27 dynamic element* of the TPSO engine, shown in figure II, coneiet of 10 
elements in the compression system and 17 elements for the remaining components of the 
engine. On* modification was performed to the configuration of NAS3-20292 involving the 
interface location of the fan elements. The boundary location of the fan diecharf* elements 
was relocated from aft of Stator 3 to forward of 8tator 3 for both the fan bypass and cor* 
streams. This allows the dynamic element boundaries to coincide with the test data probe (Sts. 
2.3 and Sts. 2.3F) locations. This dynamic element configuration with four fan elements (two 
in the core stream and two in the bypass stream) made up the baeelin* TF30 configuration. 
The model wee operated with only one-dimensional effects to enable comparison with the 
two-dimensional fan model. As will be discussed later, the four-element fan was also con- 
figured with six and eight elements to investigate fan element configuration influence on 
model response. 


Frequency response of the fan discharge pressure appears on bod is plots (gain and phase) 
in figures 12 and 13 for the fan cor* stream and fan bypass stream, respectively. The engine 
plots originate from test frequency data and the model plot* com* from the linearized transfer 
function method described in Appendix A. The baseline two-dimensional fan model response 
shows improvement over the one-diraensinal fan model. The fan core stream of the 
two dimensional modal, as shown in figure (2, achieves resonant responses, whsreas the 
one-dimensional model produced only a gradual roll-off in gain and phase. At the low 
frequencies (less than 2ft Hz), the two-dimensional model is more attenuated than engine data, 
hut achieving resonant frequencies is much more important for this extended frequency model. 
These resonant frequencies of the model at 33 and 70 Hz compare with engine data jf 30 and 
82 Hz. Kven though the one-dimenaional fan showed resonant frequencies in the bypass 
stream, shown in figure 13, the overall trend of the two-dimensional model more closely 
duplicates engine data and the large attenuation at 7 Hz has disappeared. 


Discrete Frequency Com peri eon 

The model was exercised with discrete sinusoidal frequencies of 13. 30. 60, 82, and 100 Hz 
on P . Them* frequencies were aelected because of resonant peaks and valleys in engine data. 
This study generated data for comparison of model response with large amplitudes (identical 
to engine data) to the linearized transfer function method wherein the partials were generated 
with only (U‘. steps, as shown in Appendix A. In addition to determining the non linear 
effects due to the large amplitudes, the results could be compared tit discrete frequency engine 
data. 


Surprisingly, the model results of the discrete frequency analysis closely approximated 
the linearized transfer function analysis, as shown in figures 14 and lft, with only differences in 
phase denoted for high frequencies. Only three cycles were run at each frequency, but by 
examining the transient pressure traces shown in figures 16 through 20 the startup inac 
curacies appear to be gone. Tables 1 and 2 present the detailed results of the discrete 
frequency analysis. 
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Figure 12. TF30 Engine Frequency Response of Fan Core (S»-a. 

Inlet (Sta. 2.0) for One- and Two-Dimensional Models 
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Figure 13. TF30 Engine Frequency Pesponse of Fan Bypass (Sta. 2.3F) 
to Inlet (Sta. 2.0) for One - and Two-Dimensional Models 











Figure 17 • TF30 Simulation Fan Response 
Variation (30 Hz, 8? Amplitude) 
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Figure 2(K TF3H Simulation Fan Response to Inlet Pressure Cyclic 
Variation (1H0 A* Amplitude) 
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Stability Limit* 

The final evaluation of the improvement in dynamic response due to the radially coupled 
fan involved assessment of the . lability limits of the fan. The evaluation wns performed with a 
15 H/ sinusoidal frequency on engine inlet pressure to determine the amplitude of pressure 
required to stall the engine. The results showed a dramatic improvement over the 
one dimensional model, as detailed in figure 21. Core stream stall in the simulation occurred at 
an inlet pressure amplitude of 25% (peak to peak) compared to 28.5% for engine data, 
whereas the one dimensional fan model required a t>0% amplitude to stall the engine. Even 
then, the stall occurred in the bypass stream. 

For this evaluation, pressure amplitude increased with transient simulated time, as shown 
in figure 22. This allowed phasing of the simulation to adjust before the large amplitudes were 
imposed to drive the turbomachinery pressure ratio over the stall limits. The gradual building 
ol amplitude approximated engine testing. A slower rate might have produced an even closer 
correlation between model and engine data. The simulation produced a stall at the minimum 
inlet pressure, which duplicates the engine results reported in Reference 2. 
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Instability i hi urred in ths fin rare and kw pnour* carnpraaaor whin the combined 
pressure ratio and flow were driven beyond rig still limits, u shown in figure ‘23. Figure 24 
shows that tile fin bypass stream is not near its rig still line, although this wis the stilling 
component in the one dimensional model evaluation (Reference I). 

Exceeding the stability limit of a component can be viewed as the transient pressure ratio 
exceeding some value for a finite amount of time. The transient pressure ratio during a 
sinusoidal disturlwmce can lie analytically related to the discharge pressure response (gain and 
phase! to inlet pressure. Shown on figure 2ft -A are two examples of transient pressure ratio 
excursion* as functions of exit pressure response. One example has 0 db gain, but 90* phase 
shilt. Pressure ratio increase, which is proportional to percent exit pressure minus percent 
inlet pressure, reached a maximum value of 141%. 

I’he second example has zero phase shift, but a -6.0 db gain and resulted in a 50% 
maximum increase in pressure ratio. Presented on figure 2ft B is this pressure ratio rela- 
tionship over a range of gain and phase where the increase in pressure ratio is normalized with 
inlet pressure • amplitude. The frequency response values from the engine data. 2 D model, 
and I I) model are also presented on figure 25-B for comparison. 

The two dimensional fan model very closely approximates engine data for both the fan 
hypo** and fan core/low pressure compressor streams. The lower attenuation encountered in 
thi one dimensional unwlel for the core stream, even though phase is dose to engine data, 
results in lower pressure ratio increases, which prevents the core from reaching the stability 
limit 


Number of Fan Elements 

Investigations included the dynamic configuration of the fan lor further improvement m 
frequency res|Huise. The four element baseline configuration consisted of two elements in the 
core stream und two elements in the bypass stream with the interface located between St»tor 1 
and Rotor M. As shown in figure 26, two more cc nfigurations were evaluated consisting of six 
and eight elements, respectively. Both configurations were stage -by stage representations with 
the eight dement configuration possessing different pairings of rows and a separate element 
for the space between Rotor 6 and the splitter. 

Figures 2? and 28 detail the frequency response results of the trades on element 
configuration. The tan discharge pressure response of the core stream iSta. 2. Ml and the bypass 
streum (Sta. 2.MK> appear on the same plot for e^ch of the configurations with a reference plot 
of engine test data The pressure gain as shown in figure 27. of the engine data shows much 
less attenuation in the core stream for the low frequency (less then MO Hz) region than the 
bypass stream. The baseline four element configuration produced about the same attenuation 
for both streams. The six element configuration shows separation of the two streams with the 
bypass stream exactly matching engine data. The eight element configuration shows further 
separation, but because the bypuss stream indicated too much attenuation. The resonant 
f r *quencies of MO and 82 H/ for the engine and MM and 70 Hz for the model remain about the 
same lor the different fan configurations, except that the eight element fan configuration 
7o Hz peak is very attenuut d. 

Similar trends existed for phase response between the different tan model configurations, 
as shown in figure 28. Again, the six element configuration more closely approximated engine 
data than the other configurations 
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Two Theoretical Examples 



Figure 25. Transient P/P as a Function of Exit Pressure Response 
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Figure 26. Fan Dynamic Element Breakup Configurations Investigated 
(TF30 Engine Simulation) 
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^ure 27. Fan Resnonse (Gain) Changes with Dynamic Elen 
(TF30 Engine Simulation) 












The dynamic element configuration trade study reached the conclusion that the 
six-element stage -by stage fan configuration is best suited for this application. It matched 
engine data except for the fan core stream in the low-frequency region, and in the magnitude 
of the high-frequency response amplitude peak. The low-frequency differences could be 
explained by the low pressure and/or high pressure compressor representations. 


F100 ENGINE (REMOTE SPLITTER PAN) MODEL 


The F 100(H) is an augment* 1, mixed-flow twin-spool turbofan similar in the arrangement 
of components to the TF30 engine, as shown in figure 29, Compression system differences from 
the TF30 engine include a three-stage fan without a low-pressure compressor, and a ten-stage 
high-pressure compressor. Variable inlet guide vanes are on the fan and the first three stages 
of the high-pressure compressor. The physical splitter between the core and the duct bypass 
streams is located 14 cm (5.5 in.) aft of the last fan rotor (remote splitter configuration). 
Another version of the F100(3) engine locates the physical splitter 1.3 cm (0.5 in.) aft of Rotor 
3 (proximate splitter configuration). 

Element Configuration 

The dynamic element configuration of the FIDO engine model duplicates the TF30 engine 
model, except for the compression system. The high-pressure compressor was divided into five 
elements, which represents the same number utilized for the combined low- and high-pressure 
compressors in the TF30 model. The fan was configured with eight elements, four in each 
stream as shown in figure 30. This fan configuration duplicated the best staged configuration 
of the TK30 model (six elements) with an additional element in each stream forward of the 
splitter to account for the increased axial length of 14‘cm (5.5 in.). 

Fan Characteristics 

The FIDO remote splitter fan model was matched at the fan steady-state operating 
condition for u (>. ( .) Mach number and 50, <HK) ft altitude (0.9/50k) flight condition. The dividing 
streamline purumeters (AKKA and r, as described in Appendix B) were fixed at a bypass ratio 
condition equal to 0.07. This configuration was operated transiently by slowly increasing and 
decreasing bypass exhaust pressure in order to verify fan steady-state performance character- 
istics during the dynamic mode. These very slow transients resulted in essentially steady-state 
performance. As noted in figure HI, the fan warpage characteristic matched exactly with rig 
and engine data. 

The resulting model radial and axial physical (lows after these slow transients to higher 
and lower bypass ratios appear in figure H2. The trend of radial flow in the fan rows indicates 
rapidly increasing physical flows progressing through the fan. but the corrected flow variation 
is much less. At first observation, the direction of radial flows in the turbomachinery rows 
appears to be in the wrong direction. The radial flows for the higher bypass case progress 
toward the ID. However, this direction is correct in order to represent the movement of the 
dividing streamlines, as shown schematically in figure 33. The OI) row performance character- 
istics represent a function of corrected flow per unit area (</>/A). With the area fixed from the 
initial steady state splitting streamline location, the OD stream diverts flow in the radial 
direction to obtain the (</>/A) required to produce the correct pressure rise for this increased 
bypass ratio condition of figure 33. This diverted radial (low represents the movement of the 
dividing streamline toward the II). Downstream of Rotor 3. the radial pressure warpage 
generated by the row characteristics produces a radial flow toward the OD, which represents 
the direction of the bending of the splitting streamline in relation to the initial steady-state 
splitting streamline location. 

Therefore, the following observations were made concerning this radial flow model: 

• The direction of the radial flows in the turbomachinery rows is the 
direction that the splitting streamline moves. 

• The direction of the radial flows in front of the splitter is the direction of 
streamline bending relative to the steady state location. 
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Physical Flows - Ib/sec 



(A) Zer o Cross Flow (Base Condition) - Bypass Ratio « 0.67 



(B) Bypass Duct Exit Pressure Reduced * Bypass Ratio 0.717 



(C) Bypass Duct Exit Pressure Increased - Bypass Ratio 0.63 

FD 2UJ213A 


Figure 32. PWsicnl Flow*? through Fnn ns Bypass ^uct Fxit Pressure is 
Fhan>iod (Finn Homot* 1 





Figure 33. n ynnnic Radial Flows Represent Movement of dividing 
Streamline Between Fan op and TP Stream - Bypass Ratio - 0.717 









The fan representation with radiai flow* matched steady-state characteristics (figure 31) 
when the derivatives of the dynamics were very small tslow transient) This points out the 
potential of using radial flow representations in steady-state modeling to obtain accurate 
bypass ratio effects. 


High Frequency Comparison 

This portion of the investigation compared the FIDO engine model to high frequency 
rumble tAo Hr) transient data generated on engine FX219-13 at the Arnold Engineering 
Development Center (AKIM’) in 1974. Available high-response instrumentation data tracea 
included augmentor pressure, fan exit OD and ID pressures, and main-burner presaures. 
Engine model evaluation proceeded by creating a model interface at the augmentor location to 
accept the rumhlr pressure oscillations of engine data as model input. 

The resulting model response in terms of fan discharge pressures and mainburner 
pressure to the input augmentor pressure rumble of f>2 Hz at a flight condition of 0.9/ftOk 
appears m figure 34 The fan discharge pressures of the model, particularly the OI) stream, 
show good correlation in phase, hut amplitudes are too large indicating less attenuation of the 
input augmentor pressure. 

This study reached the conclusion that the lower attenuation of the fan discharge 
pressures resulted from an inadequate bypass chief model between the lan und augmentor. The 
engine model duct representation consisted of four axial elements with a constant area, 
whereas a large cross sectional area variation exists in the real engine down the length of the 
ducts. 

When the Floo fan model was first assembled, it operated in a compression system "rig" 
mode consisting of Fl'Hi fan and compressor row characteristic* and a duct module identical to 
the TF:mi engine duct to verify model operation with augmentor pressure inputs. Operation of 
this rig model with identical augmentor rumble oscillations -suited in fan discharge pressures 
shown in tigure‘3Y As can be obsc ved. the amplitude of model fan discharge pressures closely 
approximated engine data, although slightly more attenuated, while an approximately . r > r >° 
additional phase angle existed It was anticipated that a complete FHM) engine model with a 
more representative FIDO duct (length, pressure losses, and average aren) would correct thiR 
phase angle error. As noted in the engine model version depicted in figure 34. the phase 
improved, but at the expense of n large change in gain. Other differences between the rig and 
engine models included fan exit guide vanes and an increased number of fan dynamic elements 
from four in the rig to eight in the engine model. However, the duct representation change 
produced the largest influence. An adequate duct model for this frequency band width should 
include a variable area with more than four dynamic elements. 

Model compression system parameters apjieur in figures ;iti through 41 during the engine 
simulated n2 Hz rumble transient. The large rjH'l ) variation in fan bypass ratio at the s;*littet 
reflects the excursions of the fan OD flow stream relative to the ID flow stream, as shown in 
figure 3ti. Excursions of this magnitude at this high frequency point out the need for 
aerodynamic coupling of the two streams. 

Figures ;t7. 38. and 351 sh< the transient on fan maps for mass averaged. OD. und ID 
pressure ratios, respectively. The average pressure ratio moved toward the stall line with a 
maximum pressure ratio increase of 9.2*7 and a 2*7 reduction in flow. OD pressure ratio 
increased by a maximum of 9.9 *7 . while the II) pressure ratio increased by 8.9*7 . The transient 
of the high pressure compressor basically paralleled the stall line with an approximate 
* 1.2 .7 *7 flow excursion, as shown in figure it'. 
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Figure 36. F100 Remote Splitter Fan Bypass Ratio Transient Paring A/B 
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Figure 38. FIDO Renote Splitter Fan Or Pressure Ratio Transient 
During A/ p Rumble 
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Figure 41. F100 Remote Splitter Fan Warpage During A/P Rumble 



Figure 41 illustrates the model fan pressure transient warpage characteristics along with 
the steady-state characteristic. As shown previously, steady-state model performance matched 
engine data over a much wider excursion than this rumble transient. Therefore, model 
response was different with the two-dimensional fan than if fan performance during the 
transient had been forced to match this warpage (bypass ratio) correlation. 

This investigation also compared the model with the same engine (FX219-13) during an 
additional augmentor rumble transient. A flight condition of 1.6/64K was selected because the 
fan operates at a lower operating point than at the 0.9/fi0K flight condition. The comparison 
between model fan discharge pressures and engine data appears in figure 42. As in the other 
remote splitter rumble comparison, model fan discharge pressures were less attenuated than 
engine data. 


Stability Limits 

A final evaluation induced low-frequency sinusoidal inlet pressure waves into the FlOO 
engine with remote splitter model to determine the model stability limit. At a simulated flight 
condition of 0.7/30K, 8-Hz sine waves on P, were input at four different amplitudes. As would 
t>e expected, the larger input amplitudes produced larger transient fan pressure ratios until the 
stability limit was obtained, as shown in figure 43. The condition with an amplitude of 
• 21.3% (42.5% peak to peak) resulted in model stall very close to the rig fan stall line 
detailed in figure 44. 

The frequency response (gain and phase) listed below represent fan discharge pressure 
responses to the 8-Hz discrete frequency inlet pressure disturbances. 

(lain Phase 

Fun ID 0.72db -20° 

Kan 01) t).42db 13° 

The combined effect of gain and phase on transient pressure ratio increase appears in the 
analytical relationship of figure 43 und shows good correlation with the model maximum 
pressure ratios seen in figure 43 before the model stall occurs. 
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F100 KNOMK (PAOXIMATK SPUTTtfl FAN) MOOCL 
Fan ChAraetertettM 


Ths FIDO fan modal underwent two modification* to simulaU the engine with a 
proximate iplitter. The physical splitter extended forward to within 0.6 in. of Rotor 3. First, 
the “apace” element between Rotor 3 and the aplitter wa* reduced from an axial length of 6.6 
to 0.5 in. Second, a pressure calculation due to expansion/ contraction between Rotor 3 and the 
splitter was added. The area ratios for this exDansion/contraction represent the difference 
between dividing streamline location on Rotor 3 and the nose of the splitter. Thi* calculation 
procedure wa* included in the correlation of the empirical characteristics of the TF30 
proximate splitter fan, but was not required to correlate the F100 remote splitter fan data. 

The FIDO proximate splitter fan data exhibits a steeper warpage characteristic than the 
remote splitter, as indicated in figure 46. The added area ratio calculation in the model moved 
the steady state level of warpage to match the proximate data level. The smaller space 
between Rotor 3 and the splitter provides a stiffer resistance that accounts for the steeper 
warpage slope, but the verification of this characteristic after a slow transient (as done for the 
remote splitter) was not performed due to schedule and resource limitations 

High Frequency Comparison 

Thi* evaluation compared the FIDO engine mode! with simulated proximate splitter to 
rumble data from engine 1*9*49 with close coupled instrumentation. This instrumentation 
utilized Statham transducers mounted on the engine (“dose") but some distance from the 
prolie. The lower frequency response of the Statham transducer (compared to a kulite for 
high-responsc instrumentation) coupled with the dynamics of the "close" connecting tube 
resulted in some attenuation and phase shift in the recorded rumble signals. Rumble data for a 
proximate splitter Klttti with high response instrumentation was not available. A comparison 
of the 57 Hz rumble duta of this engine to model response at a flight condition of (l.4/28k 
appears in figure 4 s. The comparison appears much d<»er than for the remote splitter engine, 
hut because of the dose-coupled instrumentation which attenuated actual amplitude ex- 
cursion*. no valid conclusions van he made between the proximate and remote splitter models 

The compression system transient excursions shown in figures 46 through 53 are larger 
than the remote splitter rumble, even though the auvmentor rumble amplitudes for each were 
both about 25% peak to peak. The fan bypass ratio transient shown in figure 49 is 32% peak 
to peak, but includes a 19% change in the reduced bypass direction The slow transient 
performance characteristics of the monel were not verified. The Ian map transients of pressure 
ratio versus flow shown in figures 51, 5°, and 53 show the stall line being approached due to a 
3% decrease in flow with only a 1.5% flow increase. This trend results from the steep (choked) 
pumping characteristic of the fan. The transient in the high pressure i» ipressor shown in 
figure 47 basically parallels the stall line. 
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SUMMARY OF RESULTS 


Thi* program created a fan model using two-dimensional modeling equations to add 
dynamic radial coupling between the core stream and the bypass stream of the fan. When 
incorporated into a twin-spool, mixed flow turbofan engine simulation, this model improved 
frequency response up to 100 H*. This modeling technology improvement was verified through 
comparison with NASA TF.10 engine high-frequency data. The new fan model also matched 
engine stability limits at low frequency (15 Hr), whereas the one-dimensional fan mode! 
required twice the inlet pressure amplitude to stall the simulation. 

With this verification of the improved fan model, a high-frequency K100 engine sim- 
ulation was constructed using row-by-row compression system characteristics. In addition to 
the baseline FlOOfil) remote splitter fan, the model fan characteristics underwent modification 
to simulate the FIDO proximate splitter fan. Upon comparison, the F100 models showed good 
correlation with high-frequency (50 to 60 Hz) rumble data. The F100 remote splitter model fan 
stability limit was verified by imposing inlet planar pressure variations at a low frequency (8 
Hz). 


Specific highlights of these results include: 

TF:i(i Simulation: 

• The improved baseline two-dimensional fan model in the engine simulation 
compared to the one -dimensional fan model and NASA engine data for 
frequencies up to 100 Hz appear in the frequency response bodie plots in 
figures 11 and 10. The new model did show resonant frequencies in lioth 
the fan core and bypass streams, even though the absolute magnitudes of 
these resonant freijuencues were in slight error. 

• The baseline four element dynamic fan could he improved slightly by use 
of a st age -by stage six element configuration as shown in figures 26. 27, 
and 28. 

• Operating the simulation with discrete, cyclic inlet pressure variations of 
the same amplitude as engine data resulted in close agreement in response 
to the more cost effective transfer function method of model frequency 
response, as shown in figures M and 15. 

• The stability limit of the new fan model agreed well with engine data at 15 
Hz, as shown in figure 21. (’ore stream stall oa.irred in the simulation at 
an inlet pressure amplitude of 22*7 peak to peak compared to 28*7- for 
engine data, whereas the old model with the one dimensional fan required 
a 60 r r amplitude to stall, and the stall occurred incorrectly in the fan 
bypass stream. 

• An analytical relationship of maximum fan pressure ratio increases due to 
variations in gain and phase of fan discharge pressure was generated. This 
relationship can lie used to quantify the combined effects of discharge 
pressure gain and phase on stability (maximum pressure ratio) limits of 
turhomachinery. as shown in figure 25 for the stability limits of the TF80 
models. 
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Fh-u Hi iunh Splitter Fan Simu'ation 


• The two dimensional fan model with radial (Iowa matched engine data 
characteristics when allowed to reach ateadystate operation after a tran- 
sient. as shown in figure 46. 

• Augmentor pressure variations (rumble) at SO Hz were imposed on the 
model and resulting Ian discharge pressures compared with engine data at 
0.9/SOk and l.ti/tMk flight conditions, as shown in figures 1)4 and 42, 
rci»|HH-tively. Phase compared well and the higher gain of the model 
proliahly resulted from an inadequate duct model. 

• Transient Ian bypass ratio variations during a .*»2 Hz rumble transient 
appear in figure Ui. The large magnitude t2H‘ , V) variation points out the 
need for dynamic coupling between the two fan streams in high-frequency 
simulations. 

• Im reusing amplitude on engine inlet pressure (I* I cyclic variations re 
suited in giuKl agreement between the model fan stability limit and the fan 
rig data, as shown in figure 44. 

Fit** Proximate Splitter Fan Simulation 

• The remote fan splitter was changed to simulate the proximate fan splitter 
configuration characteristics. The change amounted to an additional ex 
pansiiin/contraction pressure calculation term downstream of Kotor 2 that 
represented an area variation Ik* tween the location of the splitting 
streamline on Kotor and the nose of the splitter. The initial steady state 
condition helore the rumble transient matched proximate splitter fan 
warpage data, as shown in figure 41). hut other off design conditions were 
not sent led A more accurate representation of the proximate fan splitter 
would require new Kotor :< characteristics that include radial llow effects 
on warpage. but tins was beyond the scope of this contract. 

• Comparison of the proximate fan splitter simulation with a? Hz rumble 

data at a flight condition of 0.4'2Xk appears in figure . r >0. Although the gain 
of the tail diM barge pressures appear much better than the remote splitter 
model, the proximate engine instrumentation was dose coupled, not 
high response instrumentation as in the remote engine. Therefore, because 
the close coupled instrumentation dipped off the maximum and minimum 
pressure amplitudes, the exact mtalel fidelity could not be determined. 


CONCLUSIONS AND MCOMMSNOATIONS 


The following conclusions and recommendations were drawn from the modeling technolo- 
gy improvements acquired under this investigation of remote and proximate fan splitter model 
simulations: 

• The two-dimensional fan model using radial, dynamic coupling between 
the fan ID and OD axial flow streams produced significant improvements 
in frequency response. 

• High-frequency (greater than 10Hz) dynamic simulations of mixed flow 
turbofan engines require dynamic coupling between fan ID and OD axial 
flow streams to approach the correct frequency response. 

• The TF30 and FiOO engine models developed under this contract exhibit 
sufficient frequency response accuracy to help plan high-frequency engine 
testing and assist in data analyses. 

• Technology improvements in duct and augmentor models should be 
pursued for turbofan engine model applications where accurate resonant 
frequencies are required. 

• Additional modeling technology might be required, such as separating the 
variable area flow path effect from the row "actuator disk” characteristic. 

• The potential of steady-state fan characteristics using radial flow for the 
bypass variation should be explored. 


APPCMCH X a 

DYNAMIC 1 EQUATIONS AND SOLUTIONS 


DyiMNdle equations 

Thin appendix describes the equations used to simulate dynamic operation of an engine 
and to calculate frequency response. The equations express the dynamics of rotor speed 
transients, unsteady compressor aerodynamics, and gas path fluid flow. 

Rotor Dynamic* 


Simulation of engine transient operation requires accounting for the effects of rotor 
dynamic* Acceleration or deceleration of the rotors occurs when turbine horsepower exceeds 
or becomes less than compressor horsepower requirements. The basic formula used to calculate 
rotor speed changes is 


dN Aq 

dt " I, 

where: 

AQ rotor torque unbalance 

I, >- rotor polor moment of inertia 

N =*■ rotor speed 


(All 


Unatoady Compraaaor Aarodynamic a 

Incorporating unsteady aerodynamic effects accounts for the deviations from 
quasi sternly lift characteristics that an airfoil exhibits in a transient flow environment. Inlet 
airflow incidence angle changes require a finite period of time to propagate to the trailing edge 
of the airfoil. This effect on airfoil performance appears in figure 54. The difference in airfoil 
lift coefficient translates into deviations in blade and vane loss characteristics for compressor 
operation, as illustrated in figure 55. The net effect on compressor operation is an effective 
modification of compressor map characteristics in a high-frequency operating environment 

The unsteady effects are computed by applying a first-order lag to the quasi-steady 
entrophy change across the rotor or stator. The quasi-steady entrophy change comes from the 
equat ion 


As - 



H fn TR, - R Cn PR, 


(A2) 


where, TR, and PR, are the total temperature and total pressure ratios defined by the rotor or 
stator characteristics. A time constant modifies the quasi-steady entrophy change that repre- 
sents a function of the rotor or stator chord length and the axial velocity at the stator or rotor 
exit, such that 


chord lengt h 

axial velocity (A3. 
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Quasi-Steady 




Tha dynamic entropy change is obtained by appling the time constant as shown below: 




_ l_ 

rs + 1 




(A4) 


Equation A2 is rewritten to solve for dynamic pressure ratio 

PR “ exp (-f-j dnTR, - {A6) 

where temperature ratio is held constant at the quasi-steady value. 

Fluid Flow Dynamic Equation* 

The general relationships that describe one-dimensional fluid flow phenomena recognize 
three basic properties mass of the fluid, its momentum, and its energy. These represent 
conservation properties in that a change to any one of them within the boundaries of a system 
must he the result of quantities passing through the boundaries. In most basic terms, then, 
these relationship* can he expressed as 


A Mass. , =- Mass „ - Mass., (A6) 

^Momentum. , = Momentum., - Momentum..,, (A7) 

AEnergy , - Energy, - Energy ,, (A8) 


where. (’V refers to a hounded system known as a control volume or finite element; "in” refers 
to quantities entering the finite element: and "out" refers to quantities leaving the finite 
element. The introduction of the independent variable time (0 allows the formulation of 
equations that describe the properties within the finite element at a particular value of time, 
such that 


m Mass Mass 


AMass At AMass. At 
At' At 


(At» 


P Momentum Momentum 


AMuuentum.At 

At 


AMomcntum At 

At (Aid) 


E Energy Energy 


AEnergy At 
At 


AKnergv _At 
At 


lAllt 


As At approaches zero, the right-hand side of these equations represent integrals with respect 
to time of the time deviations of the incoming and outgoing quantities. These may be 
expressed as follows: 


m - 

J 'i m/.'t • dt 

( A 1 2 1 

P 

1 -d’/.it • dt 

(Alol 

E 

I oE/.'t • dt 

(Alii 
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At this point, recognize that the solutions obtained roust be in terms of the parameters 
with which fluid systems are normally identified; i.e., pressure, temperature, density, and 
velocity. The first form urns presented represent time deviations of the three basic properties at 
the boundaries of the system. The formulas that follow represent quantities of the properties 
within the finite elem ent at a particular time. 


Figure 56 ;hows a representation of the geometry within a typical finite element. The 
time derivatives of he properties, either incoming or outgoing, may be represented by: 


r»m * tir 

- -„Av-W 

§1 , ?.issai. _ vam mav _ w p 

at at at at 

— - h," 

at 1 

where, 

h, - / CpdT + 0.5 pv - J CpdT T 


(A15) 

(A16) 

(A17) 


(A18) 


The mass derivative (Equation A15) is density times the area normal to the flow direction 
times the distance traveled by the fluid per time differential. This latter quantity is the 
velocity of the fluid. Time rate of change consists of two terms, the first describing the change 
in momemtum due to mass transfer and the second consisting of the summation of forces on 
the fluid in the direction of the flow. This second term represents an application of Newton’s 
Second Law to the finite element. The energy derivative, expressed in terms of total enthalpy, 
is the sum of the specific kinetic energy of the fluid and the specific potential, or thermal, 
energy contained in the fluid. The specific enthalpy is multiplied by flow to obtain the time 
rate of change in energy at the boundary. 

At this point, the time derivatives of mass, momentum, and energy have been identified 
in terms of the parameters usually associated with fluid flow. In particular, these quantities 
have been defined at the boundaries of a finite element, or control volume, assuming that the 
flowfield is one dimensional. Spatial integrals, which define the quantities of mass, momen- 
tum. and energy with the boundaries, are given by: 


m * 

%> 

1 

f pdV - f 

A 

pAdL 


p = 

f pv dV - f 

pAvdL 


E= J 

r 0.5 pv dV + 

f pCvTdV - 

f [0.5 v + CvT] pAdL 


(A19) 
(A20) 
(A21) 

A basic assumption in the implementation of the finite element technique follows that 
the parameters of elements viry in a linear fashion from one side to the other. These formula 
for the spatial variations in the general parameter are: 


p-p„ + 


(A22) 
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m 


where, L varies from 0 to L across the element. This formula is applied to variations in 
pressure, temperature, density, area, and velocity across the element. 

Substitution of these equations into the spatial integrals allows their evaluation and 
yields the following functions for the properties within the finite element: 

m - f(p, A, AL) (A23) 

p - f(p, A, v, AL) (A24) 

E - f(p, A, v, T, AL) (A25) 


j One term still remains to be defined: the sum of the forces on the fluid exerted by its 

surrounding the flow direction ( F), which appears in the equation for the time rate of 

change in momentum. Figure 56 shows that the geometry of the finit» element plays a role in 
the assessment of the forces exerted in the fluid. The boundaries through which flow passes 
I are assumed to be normal to the flow direction. The walls of the duct form the other boundary 

of the clement. Total force consists of two parts: the force exerted by the fluid outside the 
| element arid the force exerted by the walls on the fluid inside the element. The first term is 

1 simply pressure times the flow area at the boundary. In the absence of friction or any area 

i variation across the element, there are no body forces exerted by the walls in the flow 

direction. 


The sum of forces on the fluid within the finite element then becomes: 


2F - p,„A,„ - p.*A* 


(A26a) 


Up to this point, the discussion has been limited to one-dimensional flow in ducts 
without considering the changes to the properties imposed by external mechanisms; in 
particular, changes due to the forces on the fluid exerted by the turbomachinery and the rate 
of energy transfer to the fluid in the form of mechanical work, such that 


F' = force on the fluid due to mechanical components (A27) 

dW/dt * rate of mechanical work performed on the fluid by the mechanical 
components 


The spatial integrals of Equations A19 through A21 define the three equations which govern 
fluid flow within a one-dimensional finite element These equations may be written as: 


m - (2p,„A... + 2p,.„A,.„ + p„. A„„ + p A,„) AL/6 (A28) 

= J <W,„ - W J dt (A29) 

P l(3A„, + A » p,„v„, + (A , + 3AJ p...,v.„, + (A,,. + A,J (p,„v + p,.„vj] AL/12 

- J (Wv„ - W v + P,,A,„ P A .... + F) dt (A31) 


(A30) 
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K |\ (A l/< /» + p.J 20) + A ... U>J20 + p.,/30)) 4 
V (A (P /30 4 p.„/20) 4 A , (p.720 4 p..„/5)) 4 
v„ v (A (p./JO 4 P „/15) 4 A,„ ip Jib 4 p„/10))J AL/2 4 
C'v | ( 3 A + A .) p...T.„ + (A,„ 4 3A..J p.„T„, 4 
(A 4 A .,) <p T . 4 p.„,TJJ AL/12 

* f [h, W h, JW„ + dt 


(A32) 

(A33) 


In evaluating external force and energy derivative inputs to a finite element, a concept 
known as the actuator disk is employed. In this model, a series of disks of negligible length, 
across which the total temperature and pressure change as a specific function of a specific 
component characteristic, represent the turbomachinery. The effect of these changes on the 
momentum and energy flux at the incoming boundary of the finite element become part of the 
calculations by utilizing the parameters that result from the component characteristics in the 
calculation of incoming momentum and energy derivatives. This concept appears schematically 
in figure 57. The “in" station represents the inlet of the component, while conditions that 
would exist at the component exit in the steady-state condition are assumed at the “in’" 
station. Momentum and energy derivatives computed as a function of conditions at this station 
include the effects of the turbomachinery components. On this basis, the expressions for the 
time derivatives of momentum and energy presented in Equations A28 and A29 can be 
modified in the following manner: 

W >./ 4 P ,'A„ = W,„v m 4 P,„A,„ + F' (A34) 


h, ,,,'W = h, „,W,„ 4 


dW 

nr 


(A35) 


The engine simulation describes compressor characteristics on a row-by-row basis. 
Dynamically, the model utilizes finite elements to describe compressor operation. The finite 
element configuration requires that row performance characteristics be combined algebraically 
to establish the properties associated with the finite element actuator disk. 

Solution Tochidquo 

The previous section defined three ordinary differential equations (ODE) applicable to a 
single finite element within an arbitrary flowfield. This set can be applied to each of the 
elements of the system. The resulting set of equations, when solved simultaneously, describe 
the flowfield at any point in time, given the time history of the field up to that point and its 
current lioundary conditions. The simulated response of tne system consists of a series of 
solutions at the discrete values of time required by the particular problem being studied. It 
must commence with a steady-state solution which provides the required initial values of the 
properties. 

The transient solution technique applies implicit solution methods directly to the ODE. 
A simple trapezoidal numerical integration scheme is applied to the time derivatives of the 
properties, as shown in Equations A15 through A17. This results in a set of non-linear 
algebraic equations which yield solutions at liscrete values of time in an intcrativ > manner. 
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Figure 57. Finite Flement Representation Showing the Actuator Disc 
Concept 


The concept referred to above is generally known a* a relaxative method. Relaxative 
methods rely on solving a set of implicitly defined relationships in an interative manner. This 
solution method uses a set of independent variables which permit computation of both the 
right-hand and left-hand sides of the set of governing equations (Equations A28 through A33). 
A solution is obtained when, for a particular set of boundary conditions, the values of all 
independent variables satisfy the equations. 

In applying this method to the set of non-linear algebraic equations which are used to 
model the dynamic behavior of the flowfield described by an arbitary set of finite elements, a 
set of three variables describing conditions at the exit of each element must be chosen. Two of 
these must describe the state of fluid, while the third must describe its motion. Pressure and 
temperature, either static or total, are chosen for the first two, while mass flow or velocity is 
used for the third. This set of variables is sufficient to explicitly calculate both sides of the 
three equations and, hence, the three error terms associated with them. Additional pairs of 
independent variables and errors are used to implicity define boundary conditions or con- 
straints of the system that cannot be expressed explicitly. These include the flow areas in the 
turbines and at the exhaust nozzle. A sophisticated non-linear equation solver, utilizing a 
modified Newton-Raphson approach, operates on this set of independent variables and errors 
to achieve the required solutions. 

Frequency Reepon ee Calculation 

Frequency response, defined as gain and phase versus frequency of the simulation was 
evaluated by two methods: (1) transfer functions and (2) amplitude and phase analysis of 
response to discrete frequency sine waves at P, .. The capability of these analytical methods 
was consolidated into the engine simulation and could be activated via input. 

The transfer function combined a state-variable definition of the engine with the routines 
of Seidel, as reported in Reference 3. As described in Reference 4, the state-variable method 
provided a linear description of a system using first-order differential equations ir the form: 

x(t) - Ax(t) + Butt) (A36) 

where: 

x » state vector 

u * input vector 

A * system matrix 

B = input matrix 

The output system is described by: 

y(t) - Cx(t) + Du(t) (A37) 

where: 

y * output vector 

C = output matrix 

D = output/input direct coupling matrix 


The matrices in the above expressions consist of partial derivatives determined from the 
simulation of finite differences. Each state and input variable perturb individually with all 
others held constant. 
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The male variable* used for the baseline simulation include mass, momentum, and energy 
of each dynamic element, core exhaust for dynamic mixing, and fuel flow for the engine fuel 
control. The input variable is Pt,„ and the output variables are P-n, and Pn*- Formulation of 
these variables results in matrices of the following form for a 62 nd-order system: 
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Using the computer routines of Seidel furnished by NASA-LeRC and reported in 
Reference 3, matrix algebra solved Equations A36 and A37 simultaneously to define the 
transfer function y(t)/u(t). The resulting gain and phase characteristics of the simulation were 
obtained from the transfer functions. 

The second method for defining simulation frequency response utilized direct analysis of 
a transient induced by discrete frequency sine waves applied to inlet pressure. As shown in 
figure 58, direct measurement of input pressure and output pressure wave amplitude and 
phase lag were made from the time traces of the transient. 




PAN NOW CHANACTBMSTICt 


This appendix describes tha blade row performance cbaracteriatica of the F100<3) fan, aa 
modeled under this contract Thee* blade row cbaracteriatica (from Reference 5) are derived 
from empirical F 100 fan rig data plus extrapolation* baaed on cascade correlations. A static 
pressure rise characteristic simulates the stationary (stator) element while a staff- pressure 
rise characteristic plus a total temperature rise characteristic represent the rot ng (rotor) 
element. A "dividing" streamline separates the fan core (ID) flow from the f pass (OD) 
flow. A given location of the dividing streamline, which varies with bypass . 10 , defines the 
flow areas in the ID and OD streams. The assumption made here is that now per unit area 
across the span of the fan remains constant. 

Figures 59 through 68 present the fan row characteristics for the ID flow stream and the 
OD flow stream. The performance parameters of (♦) for the static pressure rise and (X) for the 
total temperature rise represent functions of the mass flow parameter (4>/k) and the speed 
parameter. These parameters are defined in the following manner for either the ID stream or 
the OD stream, where speed parameter accounts for variable guide vane effects: 



(N//9) 

* K1 + (K2)a 


(Bl) 


(B2) 


(B3 ) 


SI’EFD PARAMETER 


(BA) 


whore: 


A - stream area (ID or OD) defined by location of dividing streamline 
K1.K2 - empirical constanta relating vane angle and speed 
N ** rotor speed 

P - pressure 

r « radius to center of stream area (ID or OD) 

T * temperature 

W - physical flow rate 

- angle of inlet guide vanes 

- exit * inlet 

• pressure/standard pressure 

- temperature/standard temperature 

Subscripts 

B - base radius at bypass ratio * 1.0 

D - design conditions 

s •* static 

T -total 















Figure 62* F100 Fan Stator 1 Pressure Rise Characteristics 
























? igure 65. 




















Figure 68. F100 Fan Discharge Strut Pressure Rise Characteristics 
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